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Abstract. Recently, we have developed a unified approach to QED in arbitrary lin- 
early responding media in equilibrium — media that give rise to absorption [Phys. 
Rev. A 75, (2007) 053813]. In the present paper we show that, under appropri- 
ate conditions, the theory can be quite naturally generalized to amplifying media 
the effect of which is described within the framework of linear response theory. 
We discuss the limits of validity of the generalized theory and make contact with 
earlier quantization schemes suggested for the case of linearly and locally re- 
sponding amplifying dielectric-type media. To illustrate the theory, we present 
the electromagnetic-field correlation functions that determine the Casimir force 
in the presence of amplifying media. 

1 Introduction 

In macroscopic electrodynamics, the effect of an arbitrary linearly responding medium in equi- 
librium can be described by means of a complex conductivity tensor of the unperturbed medium 
left to its own resources. It is well known that any electromagnetic excitation in such a medium 
eventually dissipates due to unavoidable absorption in the whole frequency range. On the other 
hand, if a medium has been brought out of equilibrium, it may act, in a certain frequency 
range, as an amplifier until after some relaxation time it will have reached the equilibrium state 
again. A detailed description of the temporal evolution of the medium properties during such 
transient processes is a highly difficult problem in general. However, when — by implementation 
of a suitable energy supply mechanism — a medium is pumped in such a way that an (externally 
controlled) quasi-stationary regime is established and maintained for a sufficiently long time, 
one may describe the effect of the amplifying medium by assigning to it a conductivity tensor. 
The properties of such an amplifying medium arc thereby regarded as being fixed — an assump- 
tion which is of course only justified in a time interval within which significant changes of these 
properties can be ignored (for the purpose at hand). Clearly, the conductivity tensor that de- 
scribes the linear response of an amplifying medium can be expected to have all the properties 
known from an absorbing medium, except that it is no longer the kernel of a positive-definite 
integral operator. 

Both absorption and amplification introduce additional noise into the system — noise that 
must be carefully taken into account in quantum electrodynamics (QED). Recently, we have 
developed a unified approach to macroscopic QED in arbitrary linearly responding media de- 
scribed by a conductivity tensor for absorbing media pp. In the present paper, we extend this 
quantization scheme to allow also for amplifying media. In Sec. [21 we outline the main aspects 
of the scheme and indicate at which points modifications are necessary when a medium be- 
comes amplifying. The extension of the theory to amplifying media is then given in Sec. [3l 
where we also make contact with earlier quantization schemes for the electromagnetic field in 
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(locally responding dielectric) amplifying media [2|3j . In this context, we also address the prob- 
lem of the necessity of a modification of the Hamiltonian, which has not been considered in 
Refs. [2|3| . It is worth noting that without this modification the field operators, whose positive- 
frequency components are associated with both annihilation and creation operators in the case 
of amplifying media, would fail to satisfy Maxwell's equations. To illustrate the theory, and 
to make contact with recent work on the Casimir force in the presence of linearly amplifying 
media [4], we present the vacuum-field correlation functions that are relevant to the Casimir 
force. The example shows that, as expected, the form of the field correlation functions changes 
significantly, compared to the familiar form observed in the case of equilibrium media, when 
amplifying media are involved — a point that has not been taken into account in Ref. |4J. In 
Sec. 21 we discuss in detail the limits of validity of the theory. Finally, we give a summary and 
some concluding remarks in Sec. 



2 Sketch of the quantization scheme 

It is well-known that Maxwell's equations imply that the electric field in the frequency domain, 
E(r, oj), obeys, in the presence of matter, the (classical) equation 

oj 2 

V x V x E(r,w) - T E(r,u) = ifj, ojj( r > oj), (1) 
cr 

where j(r, oj) is the total current density, i. e., the current density that covers all the matter (on 
a chosen length scale). If j(r, oj) can be attributed entirely to a medium whose internal atomistic 
structure need not be resolved, we may assume, within the framework of linear response theory, 
the constitutive relation 

j(r )W ) = |dV(?(r,r».E(r»+j N (r, U ), (2) 

where Q(r,r',oj) is the complex, macroscopic conductivity tensor of the unperturbed medium 
in the frequency domain |5|6j . and j (r, a;) is a Langevin noise source. Equation ([2|) covers all 
the possible features of a linear medium (in particular, temporal as well as spatial dispersion). 
According to Onsager's reciprocity theorem [516] . the conductivity tensor fulfills the reciprocity 
relation Qy(r,r',u;) = Qj,-(r',r, oj), which we will adopt throughout the paper. For chosen oj, 
Q(r, r', oj) can be assumed to be the integral kernel of a reasonably well-behaved integral oper- 
ator acting on vector functions in position space. In particular, in the spatially non-dispersive 
limit, it may become a (quasi-)local integral kernel, i.e., a linear combination of 5-functions 
and their derivatives. It should be noted that the decomposition of the conductivity tensor into 
real and imaginary parts, Q(r, r', oj) = Re Q(r, r', oj) + ilm Q(r, t',oj), corresponds, due to the 
reciprocity of Q(r,r',uj), to the decomposition of the associated operator into Hermitian and 
anti-Hermitian parts, 

<r(r, r', oj) = Re Q(r, r', oj) = § [Q(v, r', oj) + Q+(r', r, oj)] , (3) 



r(r, r', oj) = Im Q(r, r', oj) = ^ [Q(r, r', oj) - Q+(r', r, oj)] . (4) 

Throughout the paper, the superscripts T and + are used to indicate transposition and Her- 
mitian conjugation with respect to tensor indices. Since Q{v,v',oj) is the temporal Fourier 
transform of a response function [5|7|8j . it is analytic in the upper complex oj half-plane, fulfills 
the Kramers-Kronig (Hilbert transform) relations, and satisfies the Schwarz reflection principle 
Q*(r,r',oj) = Q(r,r', —oj*). Since an unperturbed medium left to its own resources is neces- 
sarily an absorbing one, the operator associated with the integral kernel er(r, r', oj) is a positive 
definite operator (see, e.g., Refs. [516] ). i.e., for any (quadratically integrable) vector function 
v(r) the inequality 

J d 3 r J dVv*(r)-cr(r,r',w)-v(r') > 0, (5) 
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which is characteristic of absorbing media, must hold. 

Inserting Eq.© in Eq. |T|), we find that E(r, w) satisfies the intcgro-diffcrcntial equation 

V x V x E(r, u) - E(r, w) - i^u / dV Q(r, r', u>) • E(r', w) = wj N (r, w), (6) 

which can be regarded as being a Langevin equation. Note that j (r, w), which is a classical 
quantity yet, has vanishing mean value (as a stochastic variable) and is related to the conduc- 
tivity tensor through a fluctuation-dissipation relation. For a detailed discussion of Langevin 
equations, we refer the reader to, e.g., Refs. |5|9|10j . The solution to Eq. ([6]) can be given in 
the form 

E(r,w) = J dV G(r,r» • j N (r», (7) 

where G(r, r', u>) is the retarded Green tensor in the frequency domain. It satisfies Eq. ^ with 
the (tensorial) (5-function source, 

UJ 2 f 

V x V x G(r, s, uj) - — G(r, s, oj) - i^uj J dV Q(r, r', w) • G(r', s, oj) = I6(r ~ s) (8) 

[I, unit tensor], together with the boundary condition at infinity, and has all the attributes of a 
Fourier transformed response function just as Q(j, r', uj) has them. In particular, it is analytic in 
the upper complex uj half-plane and the Schwarz reflection principle G* (r, r', u>) = G(r, r', — uj*) 
is valid. Since Q(r,r',uj) is reciprocal, so is G(r,r',uj), G(r, r', uj) = G T (r', r, uj), and, for real 
uj, the generalized integral relation 

Hotu J d 3 s J dV G(r, s, uj) ■ cr(s, s', uj) ■ G*(s', r', uj) = Im G(r, r', uj) (9) 

holds (see Ref. pQ). For absorbing media, Eq. ((7]) is the unique solution to Eq. ©, i.e., it is 
not to be supplemented with any solutions of the source- free version of Eq. ([6]) . 

Now, one can go over to the quantum-theoretical formulation of the theory pQ. Writing 
the operators of the electric and magnetic induction fields E(r) and B(r), respectively, in a 
picture-independent manner as 



E(r) = / dwE(r,«)+H.c, (10) 
Jo 

poo 

B(r) = / dwB(r,w)+H.c, (11) 



where the respective positive-frequency parts E(r, uj) and B (r, oj) — (iu>) 1 V x E(r, uj) are 
expressed in terms of the henceforth operator-valued noise current density as [cf. Eq. (J7J)] 

E(r,w) = ifxou J d 3 r' G(r,r',uj) • j N (r',^), (12) 
B(r,^) = Mo Vx |dVG(r,r', W ).j N (r', W ), (13) 

one can prove (Appendix [X| that the well known fundamental commutation relation 

[E(r),B(r')]=iW x Itf(r-r')/eo (14) 

is satisfied if one lets 

[j N (r,uj),j\,(r',uj')]=— ( r(r,r',uj)S(uj~uj'). (15) 
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To ensure that the temporal evolution of the electromagnetic field is in accordance with 
Maxwell's equations, the operators J N (r, lo) and jjj(r,w) have to evolve, in the Heisenberg- 
picture, like ~ e~ luJt and ~ e %U3t , respectively. Hence, a Hamiltonian H as a functional of 
j N (r,w) and j|j(r, lo) needs to be introduced such that 

[j N ( v ,u),H] = fiwl N (r,w), (16) 
which constrains the Hamiltonian to take the form 

H = nJ d w yd 3 rydVj^(r,c).p(r,r',c).j N (r', W ), (17) 

apart from an irrelevant c-number contribution. From Eq. (|16[) together with Eqs. (|15p and 
(fT7|) one can see that p(r, r',u>) has to be chosen to be the kernel of the integral operator that 
is the inverse of the integral operator associated with er(r, r', lo), 

J d 3 sp(r,s,w) • er(s, r',uj) = J d 3 s cr(r, s, to) ■ p(s, r', to) = IS(r - r'). (18) 

Note that, by means of the correspondence j(eow) -1 Q(t, t',lo) <-» x( r i r 'i w )) where x{ r i r ' '> w ) 
is the (nonlocal) dielectric susceptibility tensor, the basic commutation relation (|15[1 becomes 
equivalent to the commutation relation derived from a microscopic, linear two-band model 
of dielectric material [11], which has been used to study the quantized electromagnetic field 
in spatially dispersive dielectrics |12|13j . Provided the observables of interest — including the 
Hamiltonian — can be viewed as functionals of j (r, lo) [rather than functionals of the individual 
current contributions in a possible decomposition of j (r,oj), see Ref. [I]], Eqs. (fT5|) and (JTTJ) 
can be regarded, in view of the fluctuation-dissipation theorem(s) (see, e.g., Ref. [5]), as being 
unique, and hence, as invariable fundament of the theory. 

If, within some frequency interval, <x(r, r',w) does not correspond to a positive definite in- 
tegral operator so that the integral in the inequality ((5|) becomes negative for suitably chosen 
(quadratically integrable) functions v(r), linear amplification is possible in this frequency in- 
terval. It is not difficult to prove that the equations given above can be maintained even if the 
positivity condition J^J is abandoned, provided that nevertheless (i) the Green tensor retains 
its analytic properties and (ii) p(r,r',u) continues to exist (for a detailed discussion of these 
conditions, see Sec. [4]). However, the diagonal form 



H = J duhuj d 3 rft(r,w) -f(r,w) 



(19) 



in which the Hamiltonian (|17|) can be brought in the case of absorbing media by means of a 
linear (and invertible) transformation of the variables, 

i 

j N (r, to) = (^-j 2 J dV K(r, r', «) • f (r\ lo), (20) 

cannot be maintained if the condition ^ is abandoned. In Eqs. (jl"9"]) and (|20p. f(r,oj) is a 
bosonic field, 

[f(r, w), ft (r', J)] =IS{T-r')S{io-u'), (21) 
and the integral kernel K(r, r', u>) has to obey, for real lo, the integral equation 



J d 3 s K (r, s, lo) ■ K+{r', s, lo) = cr(r, r', lo) 



(22) 



for Eq. (fl"9|) to be an equivalent representation of Eq. (|17| . In the case of absorbing media, one 
may evidently construct a possible integral kernel K (r, r', lo) in the form of 

K(t, t',lo)= J da a? (a, lo) F(a, r, lo)F* (a, r', lo) (23) 
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[er 1 / 2 ^, uj) > 0], where the complete and (5-)orthonormal functions F(a,r, uj), 

J dQF(a,r,u)F*(c«,r» = IS(r-r'), (24) 

J d 3 r F*(a, r, uj) ■ F(a', r, uj) = 8{a - a'), (25) 
are defined by the eigenvalue problem 

J dV tr(r, r', uj) ■ F(a, r', uj) = <r(a, uj)F(a, r, uj), (26) 

Here, the real uj plays the role of a parameter and a stands for the set of discrete and/or con- 
tinuous quantities needed to label the eigenfunctions. (An a-integration therefore symbolizes 
multiple summations and/or integrations over all those quantities.) Note that the kernel (|23p 
is not unique but is determined by Eq. (|22|) only up to a unitary transformation, which corre- 
sponds to the possibility of redefining the dynamical variables f (r, uj) ; for details, see Ref. pQ . 
Equation (f2"2"| is inconsistent if <x(r, r',w) is not the kernel of a positive definite integral oper- 
ator so that in this case no valid kernel K(r,r',u>) exists at all. Hence Eq. (fT9|) would fail, if 
f(r, uj) were regarded as being a bosonic field satisfying the commutation relation (f2"Tj) . 



3 Extension to linearly amplifying media 

In order to include in the quantization scheme linearly amplifying media, we first note that, 
although Eq. (|2"2"|) cannot be satisfied anymore so that Eqs. (fT!?|) - (f2"l"j) do not apply either, 
Eqs. (|24p (|26|) remain valid, with a(a,uj) not being restricted to positive values anymore. We 
may therefore expand p(r,r',uj) as 

p(r,r',uj) = J daa' 1 (a,uj)F(a,r,uj)F*(a,r',uj), (27) 

which enables us to rewrite the Hamiltonian (fl~7|) as 

H= / duj huj I dasgna(a, uj) g' (a, uj)g(a,uj), (28) 



where 

g{a,uj) = {^j 2 \a(a,uj)\-i J d 3 rF*(a,r,uj) ■ j N (r,w). (29) 

With the help of Eqs. (15]), (25]), and (TSHJ) , it is not difficult to see that 

uj),g\a', uj')] = sgn a(a,uj) S(a - a')S(uj - uj'). (30) 

The g(a,uj) may be viewed as non-bosonic generalizations of the natural variables considered 
in Ref. [Tj, as they coincide with the latter in the purely absorbing case where sgner(a,a;) = 1. 
Inversion of Eq. (|29|) by means of Eq. (f24|) yields 



l N (r,uj)=f^\ J da\a(a,uj)\s F(a,r,uj)g{a,uj). (31) 

The commutation relation ([3H)l shows that g{a,uj) [<7^(a,w)] is a bosonic annihilation (cre- 
ation) operator for positive eigenvalues a(a,u>) whereas for negative ones, g(a,uj) [g^(a, uj)] is 
a creation (annihilation) operator. It makes therefore sense to rename the operators according 
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to this behavior. Thus, denoting in each of the two cases sgna(a,oj) — ±1 the respective an- 
nihilation operator by b{a, oj) and the respective creation operator by w(a,uj), we can rewrite 
Eqs. (JUJ), dHDJ, and (j3"T|) as 



H = J du> huj J da b\a, u>)b(a, ui) — J duhuj J da b{a, uj)b^(a, uj), (32) 
(+) ' " (-) 



[b(a,uj),P{a',uj')) = 5{a - a')S(uj - to'), (33) 



and 



j N (r,w) = ^— ^ |y da <r 2 (a, F(a, r, uj)b(a, ui) + J da [— a(a, ui)] 2 F(a, r, u)P (a, uj) 

(+) (-) 

(34) 

respectively, where the notation Jy_> da ■ ■ ■ means that the integration extends over those values 
for which sgn a(a,uj) =±1. Note that the ranges of integration depend on the chosen frequency 
in general. It may be convenient to change to normal order the second term in Eq. (|32p . 
b(a, uj)¥(a, oj)^ : b(a, ui)w(a, oj) := b'(a, oj)b(a, uj), i. e., to replace the Hamiltonian in Eq. 
with 

a,Lu)b(a,uj), (35) 



H = J dwhhjj da sgn a (a, uj) b^ ( 



which differs from the Hamiltonian in Eq. (|32[) by an (infinite but) irrelevant c-number. Since 
the state space of the system is to be constructed by means of the bosonic variables b(a, uj) and 
w{a,ui), the use of Eq. ([33)) in place of Eq. (|3"2"|) is equivalent to a redefinition (renormalization) 
of the zero of energy by the condition of absence of quanta (vacuum state). It should be stressed 
that the temporal evolution of the variables b(a,uj) and b*(a,uj) which follows from Eq. (|3"5"]) 
[together with Eq. ([33]) ] is sensitive to the sign of a(a,uj) in just such a way that Eq. (f34|) 
always represents the positive-frequency part of the noise current density, as required. 

From Eq. (|35p it is seen that there is a continuum of negative energy eigenvalues in the 
case of amplification, in addition to the positive-energy continuum associated with absorption. 
Thus, the vacuum state can no longer be said to be the ground state of the system — there is 
in fact no ground state as the continuum stretches down to — oo. This somewhat unpleasant 
feature is due to the fact that the pump mechanism which prepares the medium to act, in 
some frequency interval, as an amplifier is not dynamically included in the theory. Clearly, the 
approximation to treat the effect of pumping within the framework of linear response theory 
breaks down when states with large negative energies are significantly involved. 

By means of the transformation 

= J da F (a, r,uj)b(a,uj) + J da F(a, r, uj)tf{a, uj), (36) 

(+) (-) 

vectorial field variables f (r, uj) can be introduced, which can be viewed as generalizations of 
the variables f(r,w) introduced in Eq. (j20[) for the case of purely absorbing media. Employing 
Eq. ([25]) to invert (the first equation in) Eq. ([36]) and inserting the result in Eq. ([31]) , we obtain 

the generalization of Eq. ([20]) [f(r, uj) i— > f(r, u>)], from which we can read off, by comparison 
with Eq. ([20]) . the generalization of the kernel K(r, r', uj), in the form of its eigenfunction 
expansion. Not surprisingly, it is given by Eq. (f23|) with a(a,uj) being replaced by |cr(a,w)|. 
(It is of course possible to consider equivalent kernels just as in the case of purely absorbing 
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media.) Unfortunately, the variables f(r, w) do not diagonalize the Hamiltonian in general, and 
are thus less useful than in the case of purely absorbing media. From Eqs. (f30|) and (|36|) . it 
follows that 

$(r,w),fr(r',u')]=5(u-u') J da sgncr(a, w) F(a, r, w)F*(a, r', w). (37) 

The integral appearing on the right-hand side is the kernel of a parity-type operator, which 
reduces to the unit operator in the case of purely absorbing media, as it should be. 

A noteworthy simplification occurs if the medium response is strictly local — an assumption 

that is typically made in the study of amplifying media. In this case, the f (r, u>) are related 
to the g(a,u>) [a t— > (i,r)] in a very simple way just as in the case of purely absorbing media 
that respond locally [since the eigenfunctions F(a,r, w) are spatially localized in this case, see 
Ref. pQ]. Focusing, for simplicity, on isotropic media, we may then rewrite Eqs. (f31~j) and (|28f . 
respectively, as 

j N (r,o;)=('^ 2 |a(r ! a;)|if(r, W ), (38) 

and 

H = J duhuj d 3 rsgaa(r,u)V(r,u) -f(r,w), (39) 

and ([37|) simplifies to 

[i(r,u),fr(r',uj')] =sgna(r,u)I6(r-r')6(LJ-Lj'). (40) 

Now we can switch to genuine bosonic variables by renaming f (r, u) as b(r, u>) and b^r, u>) for 
sgn<r(r,[x)) = 1 and sgncr(r,o;) = —1, respectively, so that Eq. (|4"0")) changes to 

[b(r, w), bt(r', </)] = W(r - r')5(w - a/). (41) 

The noise current density in Eq. (|34p and the Hamiltonian in Eq. (|35p . respectively, can then 
be expressed in terms of b(r, ui) and b^(r, oj) as 

i 

1,(t,l>)= r^j 2 {^^(^^Kr.^+^-^r.^h^r.w^bt^a;)} (42) 

[#(x), unit step function] and 

F = ^ dujhuj d 3 rsgaa(r,aj)b^(r,uj) ■ b(r,w). (43) 

Taking into account that tr(r, w) can be related to the imaginary part of the dielectric permit- 
tivity according to a(r,u) — Equj Ime(r, u>), we find that Eq. (|42| is nothing but the equation 
suggested in Ref. [5] for the case of isotropic, locally and linearly responding dielectrics that also 
allow for amplification. Note that, as already mentioned for the more general case of Eq. (l34l) . 
both the term associated with h(r,u>) and the term associated with br(r, a;) in Eq. (|4"2"|) give 
rise to positive-frequency parts of the noise current density. 

As an illustration, let us consider the correlation functions (E(r)E(r')) and (B(r)B(r')) for 
the case where the system is in the vacuum state |0) [b(a,w) |0) = 0]. As well known, they 
play an important role in the calculation of the Casimir force (see, e.g., Refs. |14|15j ). Using 
Eqs. ([TO]) — <pT5]) together with Eq. (fM)) . taking into account Eq. (|33[) . and recalling the reciprocity 
of the Green tensor, one can show by straightforward calculation that 

(E(r)E(r')) = ^ J°°duiui s J d 3 s J d 3 s' G(r, s, u) ■ cr av (s, s', w) • G*(s',r',w) (44) 
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and 

(B(r)B(r')) 

= J^dujujV x |d 3 S |dVG(r,s, U ).a a¥ (s,s».G*(s',rV)xV', (45) 

where the kernel <x av (r, r',oj) reads 

<x av (r,r» = J da \a(a, w)|F(a, r, lu)F* (a, r', u). (46) 

In the case of purely absorbing media, tr av (r, r', u>) is nothing but <x(r, r 1 ,oS), and by applica- 
tion of the integral relation (J9j> , Eqs. (|44|) and (j45|) reduce to the standard (zero-temperature) 
fluctuation-dissipation relations. In contrast, if amplification is allowed for, Eqs. (j44l) and (|45|) 
differ from the standard formulas by terms containing the non-vanishing difference between 
<T a v(r, r', oj) and <r(r, r',uj). In a recent study of the (zero-temperature) Casimir force in the 
presence of linearly amplifying (left-handed) material [3], the necessity of this correction has 
not been taken into account. In fact, in Ref. |4J, a (locally responding, isotropic) amplifying 
magnetodielectric has been introduced, by simply using the formulas for the correlation func- 
tions that are valid for absorbing media and replacing therein the positive imaginary parts of 
the susceptibilities by negative ones. As we have just shown, this is wrong. 



4 Range of validity 

It remains to specify in more detail the conditions that must be satisfied to apply the quan- 
tization scheme. Let us first consider the question as to whether the Green tensor G(r,r',uj) 
remains analytic in the upper complex uj half-plane if linear amplification is allowed for. In the 
case of absorbing media, it is well-known that solutions to the source-free macroscopic Maxwell 
equations [i.e., solutions to the homogeneous version of Eq. (|6|)] for real frequency uj can be 
ruled out because of their divergent spatial behavior. The same is obviously 'even more' true for 
frequencies in the upper complex uj half-plane. Since the existence of permissible solutions to 
the source-free equations is known to manifest itself mathematically in the form of singularities 
of the Green tensor, only the lower complex uj half-plane is a possible location of such singular- 
ities in the case of absorbing media. In contrast, for linearly amplifying media it may happen 
that permissible solutions to the source-free macroscopic Maxwell equations exist even if u> is 
chosen in the upper complex uj half-plane. In this case, singularities of the Green tensor in the 
upper complex uj half-plane would exist, which would invalidate the proof of the commutation 
relation in Eq. HU). 

It is not difficult to imagine that such 'unwanted' solutions to the source-free macroscopic 
Maxwell equations could arise, e.g., if waves were allowed to propagate through an infinitely 
extended region of amplification. Since such regions do not exist in practice, their necessary 
exclusion from consideration is of no practical relevance. More importantly, the same problem 
of 'infinite amplification length' can also occur in the case of a region of amplification of finite 
extension, if waves can pass through the region repeatedly (due to multiple reflections) with a 
net gain per round-trip. Therefore, setups where an amplifying medium is part of an arrange- 
ment of bodies that act as a high-Q resonator must possibly be excluded from consideration. In 
such cases, the Green tensor fails to be a causal function in the sense of linear response theory. 
Of course, what breaks down is not really the principle of causality but the very concept of 
linear amplification — as fields with higher and higher energy develop, the non-linear dynam- 
ics can no longer be disregarded. In fact, there is no need for extra criteria to exclude from 
consideration setups that would mathematically support 'unwanted' solutions — such setups are 
already excluded implicitly by the assumption that the approximate concept of linear amplifi- 
cation is applicable. Hence, by this basic assumption, the Green tensor is forced to be analytic 
in the upper complex uj half-plane, and all the other important properties of the Green tensor 
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(in particular, its high- and low-frequency behavior and its decay behavior for large difference 
of the spatial arguments) are the same as in the case of absorbing media. 

Next, let us answer the question of the existence of the inverse of the integral operator 
associated with <x(r, r', to) in the case of linear amplification, i. e., the question of the existence of 
the kernel p(r, r', to) according to Eq. (Tl8|) . Obviously, p(r, r', to) would fail to exist if eigenvalues 
of the integral operator associated with er(r, r', to) were exactly equal to zero. Since amplification 
is limited to a certain frequency range, and since frequency is a continuous variable, each 
eigenvalue a (a, to) on the negative side of the eigenvalue spectrum can be made to move to the 
positive side by tuning the frequency, so that zero eigenvalues are possible. However, it should 
be emphasized that in the case of a continuous spectrum this problem is in fact harmless. 
It is generally true that zero occurring as a continuum eigenvalue does not really preclude 
the inversion of (Hermitian) operators. A familiar example is provided by the free-particle 
Schrodinger equation considered in the whole space. There, the continuum of (5-normalizable 
(plane-wave) eigenfunctions includes the (spatially constant) zero-energy eigenfunction, but 
this does not imply the nonexistence of the inverse of the free-particle Hamiltonian but merely 
the unboundedness of the inverse operator. Although there are good reasons to believe that in 
practice the eigenvalues cr(a, to) form a continuous spectrum, it seems advisable to have also a 
method in store to handle discrete zero eigenvalues, recalling that the calculation of physical 
quantities we are dealing with generally involves space and frequency integrations. As long as 
such a quantity remains meaningful as a whole, it does not really matter if p(v,v' ,to) is not 
literally well-defined at individual frequencies. Thus, at least on the level of final physical results, 
the effect of a discrete zero eigenvalue cannot be significantly different from that of a very small 
but non-zero eigenvalue so that any method of regularizing p(r, r', to) that implements this idea 
can be expected to give the same final results. If this is true, the problem that p(r, r', to) might 
possibly fail to be literally well-defined should be irrelevant in practice. One may simply perform 
all calculations for a class of amplifying media for which there are no problems with p(r, r',to). 
Results obtained in this way, if they are physically understandable, can then be expected to 
hold also without restriction. 

Finally let us address the problem of the unboundedness from below of the energy eigenvalue 
spectrum in the case of linear amplification. Because of the lack of a lower bound, the system 
could evolve into states of lower and lower energy by the creation of quanta in the frequency 
interval where sgna(a,tu) = — 1, in which case the theory would gradually become unrealistic. 
If the system described by the Hamiltonian in Eq. (f35|) is coupled to a second system (e.g., an 
atom), another aspect of the problem is that the second system might (but need not) become 
more and more excited, which of course also becomes unrealistic at some stage. However, such 
catastrophes could only occur if the theory were used beyond its range of validity. In fact, 
they would indicate nothing but the breakdown of the concept of linear amplification. As 
long as the concept of linear amplification applies, the unboundedness from below of the energy 
eigenvalue spectrum may be regarded as being a purely formal drawback rather than a real one. 
Needless to say that this unboundedness prevents one from constructing the canonical density 
operator — the system cannot thermalize. Therefore, it should also be clear that electromagnetic- 
field correlation functions and fluctuation-dissipation relations in the familiar form known for 
absorbing media are not applicable to amplifying media, not even at zero temperature, i.e., 
to the vacuum state. Nevertheless, all the correlation functions required can be calculated 
straightforwardly for any well-defined quantum state, in particular, for the vacuum state, as 
has been demonstrated at the end of Sec. [3] for the particular correlation functions needed for 
studying the Casimir force in the presence of linearly pumped media. 



5 Summary and conclusion 

We have shown that and how the very general quantization scheme developed in Ref. [1] for 
the macroscopic electromagnetic field in arbitrary linearly responding media in equilibrium can 
be extended in a rather natural way to include media which (in some frequency interval and 
some spatial region) are weakly pumped so that the effect of pumping can be approximately 
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described within linear response theory. In this sense, we have described the medium by a 
complex conductivity tensor Q(r,r',w) and allowed for negative eigenvalues of the integral 
operator corresponding to its real part <r(r, r',u/). Therefore, the basic condition to apply the 
theory is the validity of the concept of linear amplification for the respective problem under 
study. It has turned out that when considering amplifying media more care and prudence is 
needed than in the case of robust equilibrium media, which only give rise to absorption. 

Making contact with earlier work, we have shown that in the special case of amplifying 
linear media that can be regarded as being isotropic and locally responding, the results in 
Refs. |2l3j can be recovered. To illustrate the theory, we have calculated the electromagnetic- 
field correlation functions that determine the Casimir force in the presence of linearly amplifying 
media. As expected, they are quite different from the well-known correlation functions in the 
case of absorbing media. Since the latter are wrongly used in Ref. [4] to study the effect of 
linearly pumped magnetodielectrics on the Casimir force, with special emphasis on left-handed 
metamaterials, the results presented therein are questionable. 

For a deeper understanding of the macroscopic theory considered in the present paper, 
it would certainly be advantageous to have also available a more microscopic model of the 
quantized electromagnetic field interacting with linear media that also allow for amplification, 
especially an analog of the well-known Huttner-Barnett-type harmonic-oscillator models fre- 
quently used to study the quantized field in absorbing media |16|17j . Such a model should 
include a reservoir as in the case of absorbing media, but presumably also a second, 'inverted' 
reservoir capable of being prepared in a (formal) negative-temperature state. To our knowl- 
edge, Huttner-Barnett-type models that aim to incorporate amplification have unfortunately 
not been developed. 

We acknowledge discussions with M. Fleischhauer about the Casimir force in the presence of linearly 
amplifying media. 



A Proof of Eq. (fl4h 

Using Eqs. (TTOj) (Tl3]) and recalling the reciprocity of G(r, r', oj), we may write 

[E(r),B(r')] = -1$ J dujuj J du>' J d 3 s J dV 

x {G(r,s,u)- [j N (s,cu),g(s', w')} ■ G*(s',rV) 

+ G*(r,s,w) • [jk(s,u;)J N ( S ',u/)] • G(s',r>')} x V'. (47) 

Applying the commutation relation (|15[) and employing the reality and reciprocity of <x(r, r', ui), 
we may carry out one of the frequency integrals to obtain 



[E(r),B(r') 



Re / duu 2 d 3 s / d 3 s'G(r,s,w) • tr(s,a',u) ■ G*(s',r',u) x V'. 



(48) 



By means of the integral relation ©, we can now evaluate the spatial integrals. We find, on 
recalling the relation G* (r, r',co)= G (r, r', — v*), 



[E(r),B(r') 



duj uj Im G(r, r', ui) x V 



dwwlm G(r, r',u) 



x V'. 



As, within the framework of macroscopic QED, the presence of any medium is irrelevant at 
sufficiently high frequencies, the leading asymptotic behavior of G(r, t',uj) for frequencies in 
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the upper uj half-plane (including the real axis) with large absolute values is given by 



just as for the Green tensor in free space. Together with the fact that the most singular term of 
G(t,t',u>) at uj= is — c 2 £(r, r')/w 2 , where L(r,r') is a real tensor that is purely longitudinal 
from both sides [and thus does not contribute to Im G(r,r',w) for real u>], this asymptotic be- 
havior shows that the integral in the square brackets in Eq. (|49[) converges. It may be evaluated 
by contour-integral techniques as 



Here, the principal-value (V) integral has been changed to an integral over a contour C that 
consists of an infinitely large semi-circle in the upper half-plane (traversed clockwise) , plus an 
infinitely small semi-circle (traversed counter-clockwise) that avoids the origin in the upper 
half-plane. Note that the sub-leading (weaker than lo~ 2 ) singular terms of G(r,r',u>) at u> 
= do not contribute, irrespective of the actual nature of the singularity. Inserting Eq. (j51"j) 

in Eq. (|49|) and using X(r,r') x V = 0, we arrive at the desired Eq. IT4| [V x IS(r — r') 

= -IS(t-t') x V']. 
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